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Abstract. Traditionally the exponential canonical distributions of Gibbsian statistical 
mechanics are given theoretical justification in at least four different ways: steepest desecent 
method, counting method, Khinchin's method based on the central limit theorem, and maximum 
entropy principle of Jaynes. Equally ubiquitous power-law canonical distributions are shown to 
be given similar justification by appropriately adopting these formulations. 
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This is a report on our collaborative work done in the recent past [1-5] on an inquiry into the 
microcanonical foundation for power-law canonical distributions that occur so often in many 
physical contexts, in much the same way as the exponential canonical distributions of Gibbs do. 
There are at least two ways of understanding the exponential canonical distributions based solely 
on the principle of equal a priori probability: steepest descent and counting methods. The third 
way is due to Khinchin [6] based on the central limit theorem. Finally, the fourth is maximum 
entropy principle of Jaynes [7]. The third method is perhaps the most satisfactory 
mathematically, because only probabilistic principles are employed in establishing statistical 
mechanics of systems. In some way, the latter two methods have also the principle of equal a 
priori probability implied in their foundation. 

We elaborate on the procedures needed to deduce the power-law canonical distributions. To 
examine the possibility of obtaining the power-law-type distribution, the form of the function we 
consider is 

e q (x) = [l + (l- q )xf- q) , (A) 

which will be referred to as the ^-exponential function. This is the inverse function of the q- 
logarithmic function 

ln,(*) = n (B) 

l-q 

Throughout this paper, q is assumed to satisfy the condition 

q>\. (C) 

In the limit q—>l, e q {x) and ln ? (^:) converge to the ordinary exponential and logarithmic 

functions, respectively. For large x, the ^-exponential function asymptotically exhibits the 
power-law behavior 

2 



Our aim here is to obtain this form for the canonical distribution function by various methods of 
statistical mechanics mentioned above. 



Steepest Descent Method 

Consider a classical system s and take its N replicas s v s 2 ,L ,s N . The collection 
s = 2 l n 1S re f erre( l to as a supersystem. Let A a be a physical quantity associated with 

the system s a . This is a statistical random variable and its value is denoted by a[m a ) , where m a 
labels the allowed configurations of s a . The quantity of interest is the average of {^ a } l2L N 

over the supersystem: (1/^)^ t \ ■ ^ n microcanonical ensemble theory, the probabilities of 

finding S in the configurations in which the values of the average quantity lies around a certain 
value a [8]: 

\M\<£, (1) 



where 



M = -^£a(m a )-a, (2) 



£~0(N~ 1 - 5 ) (S>0). (3) 



The condition in eq. (3) will be revisited later. We restate eq. (1) as the step function, 
0(x) = 1 (x > 0), (x < 0). The microcanonical probability P(m 1 ,m 2 ,...m N ) associated with eq. 

(1) is 
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P(m v m 2 ,...m N )°c 9(e-\M\). (4) 

Let us select s l as the objective system. The probability of finding it in the configuration m, = m 
is given by 

p(m)= ^ P(m,ffl 2 ,...m ll ). (5) 

m 2 ,...m N 

We can represent the step function by 



B( X )= J df-f^, (6) 



as long as ji satisfies 

1 >0. (7) 

Equation (6) reduces to the traditional representation when q—>l. The condition in eq. (7) will 
be discussed later in the context of the steepest descent approximation. 
Now, using eq. (6), we have in the leading order in N 



^ + r°° e (dye) N ( 1 A 



In obtaining this equation, the q -exponential function is factorized as e q {<p\e - M]) ~ 
e q {<p e)e q {-<pM) and then e q {-<pM) is further factorized into the product over the supersystem. 
Since both \M\ and e are of 0(N~ l - s ) with 8 > 0, such a manipulation will be justified in the 

subsequent discussion of steepest descent method in the large- /V limit (as well as in the section 
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on the generalized central limit theorem). Using 0(s - \M\) = 0(s - M) - 0(-S - M) and 
performing a change of the integration variable as <fi — > Nip , we have 

£•_,•«, ^ a=l 

where sinh 9 (;t) = e ? (-x)]/2 and 



jS'=£. (10) 
TV 



Let us examine the condition in eq. (7). It is rewritten as 

l-(4-l)/rAf|±£-M| max >0. (11) 

This is a rectangular distribution function with a very narrow support. In fact, | ± £ - M| max is a 
quantity of o(N~ l ~ s ) with 8>0. Therefore, f can be an arbitrary positive constant in the 

large- N limit. 

In the leading order in N, the probability is 



P q ( m ) = p ( m > m 2> L > m A<) 



m 2 ,L ,m N 



1 P T sinh (M/>£) , r , x ,v 

x l iRHhK)-*]) 



m 2 ,L ,m w a-2 

I sinh (Atye) e (-<j)\a(m)-a]) r , ,, 

= 4 I d 4> / t „ J; exp iVlnZ g I 



(12) 



where W is the number of possible configurations allowed by the condition in eq. (1) and 
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m 



(13) 



Also, W is given by 



W= j# ^-^exp[iVlnZ 9 (0)], (14) 

ji -loo T 



because of the normalization of p (m). Using the real part /T of <fi, the steepest descent 
condition yields, with the superscripts "(5)" to indicate steepest descent method, 



dZ (S) 

— *--<), (15) 



which gives rise to 



p?w = z^ e <(-PW ~ "D ' (16) 

fl=^P,(m)a(m), (17) 



simultaneously, where P,(m) is given by 



In the field of thermodynamics of chaotic systems [9], P (m) in the above equation is referred to 

as the escort distribution. It is another probability distribution associated with the original 
distribution p q {m). These results follow from eq. (A) and de q (x)/dx = [e ? (x)J ? . Clearly, all the 

discussions become the ordinary ones in Gibbsian canonical theory with the exponential 
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distribution function in the limit g— >l + 0. The distribution function in eq. (16) is seen to 
behave asymptotically as 



p q (m) ~ 



1 



[a{m)\ 



(19) 



which exhibits the desired power law. 

The steepest descent condition leads to the fact that the arithmetic mean of {^( m a)} a _ 12L N 

coincides with the generalized expectation value with respect to the escort distribution as in eq. 
(17). 

It is worth pointing out the appearance of the Legendre transform structure [10] in this 
development. Let us introduce 



(20) 



The maximum of Z ((j)) in eq. (13) is found at the stationary point along the steepest descent 
path: 



d(3* 



= 



dp 



-a 



(21) 



In addition, we also have the equivalent result 



(22) 



provided that a = (A(m)) = P q ( m ) a(m). 
Next, we examine counting method. 
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Counting Method 

The number W(m) of configurations characterized by eq.(l) where we now isolate the mth 
state, is written as 



Ua{m)-a] + ^- £ [a(m )-a] 

^ V - /V a = 2,h ,N 



< £. 



(23) 



An important point is that the evaluation of W(m) is context dependent, in general. For example, 
it is known [1, 11] that the counting rule in the fractal space is not unique. In the present context, 
we explore the possibility of using the ^-logarithm of the density of configurations satisfying eq. 
(23). Thus, the number of configurations, Y N (a)± obeying the condition 



^ L [«K)-«] 

ly a=2,L ,N 



< £ 



(24) 



behaves according to 



\n q Y N (a)^S^(a) 



(25) 



where S ( q C) (a) will be identified below with the entropy in accordance with the definition of 
temperature which will be denoted by f5 . Here, the superscripts "(C)" are used to indicate 
counting method. 

When jV is large, the required number of configurations, W(m), is now expressed in terms of 
Y N (a) in the form: 



\n q W(m) = \n q Y N [a - ^[a(m)- afj 
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= \{Y N (a)) - [a(m) - a] « " . (26) 

da 



Let us define 



P' = [YAa)] q - l p = [YAa)] q - ldS j^. (27) 

da 



Then, equation (26) , upon using the relation in eq.(A), may be expressed as 

* = t! (-0'W-«|. (28) 
Since the total number of configurations is given by W = ^ W(m), we have, upon defining 



(/}) = W/Y„(a) = £ e< (-/3'[a(m)-a]), (29) 



Now, we turn our attention to the development based on the generalized central limit 
theorem. 



Method Based on Generalized Central Limit Theorem 

In the preceeding two sections, we considered an arbitrary discrete physical quantiy A a . Here, 

we treat a continuous quantity, the energy of the system as a representative, in order to make use 

of the limit theorems in probability theory of continuous variables. Given a system composed of 
a large number of subsystems with energies {e ; = EJB N > 0} 2 l , where B N is a positive TV- 
dependent factor to be determined subsequently, let us consider its total energy 
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E = £ x + £ 2 +L +£ N . Recall the ordinary central limit theorem. If each of E t 's obeys a common 
distribution /(F,) with ordinary finite second moment, (F, 2 ^ = | dE i F ( . 2 /(F ; ), then its Af-fold 

convolution, B N (f*f*L *f)(B N E), approaches the Gaussian distribution in the limit of large N, 
where {f*g)(x) = j dx f(x - x') g(x'). This property was exploited by Khinchin [6] to 

establish Gibbsian canonical ensemble theory within the framework of probability theory. 

Here, we note that the power-law distributions have no finite ordinary moments. However, the 

generalized central limit theorem states that the Af-fold convolution of such a distribution still 

converges to a limit distribution. The associated limit distribution is the Levy-stable distribution, 
denoted here by F (F,.). Since the energies are bounded from below, it suffices to consider the 

Levy distribution in the half space, whose characteristic function is known to be of the form [12] 



X 7 (t) = ]dE t e iE ^F Y {E) = exp -a\t\ r exp 



' sgn(0 



On 



(31) 



where a is a positive constant, y e (0, l) is the Levy index, 9 a constant satisfying |0| < y, and 
sgn(/) = //|/| the sign function of t. If each of the E ( 's obeys F (F ; ), then their A^-fold 
convolution B N [F y * F y *L *F y }(B n E) has the same characteristic function as that of the original 
F r (F,), that is, x r N \t) = [Z r {t/B N )] N = X r (*)> provided that B N is chosen to be 



B N = N l ' r . 



(32) 



F r (F ; ) is found to exhibit the following power-law behavior for large values of E-. 



(33) 



As an example of the distribution without finite ordinary moments, let us consider a power- 
law distribution 
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f(E,)=^[m+E,r (i<*<2), 



(34) 



s-1 



(35) 



In the above, is a positive constant depending on a parameter The range of s is 

determined by requiring that /(£, ) is normalizable and all of its ordinary moments be divergent. 

The logarithm of the characteristic function of the iV-fold convolution of this distribution is 
found to be [4] 



In XV) = -a \t\ exp 



isgn(f) 



On 



(36) 



Here, L(a), a, and 9 are given by 



L («) = |^( e " y - 1 )< (0<«<1), 

o y 



(37) 



a = (5-l)^ s - 1 (jS)|L(5-l)|>0, 



(38) 



e = i-s, 



(39) 



respectively. M w is a quantity, which vanishes in the limit N — > <». Comparison of eqs. (31) and 
(36) leads to the identification of the Levy index to be 



7 = 5-1. 



(40) 



To fully characterize the limiting distribution F '(E), it is necessary to relate the constant a 
to a certain statistical physical quantity. In contrast to Khinchin's discussion, we have no finite 
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ordinary moments of the distribution. This forces us to seek a different approach to determine the 

constant a appearing in eq. (36). 

In spite of the fact that the ordinary first moment, (E t ) , is divergent, the generalized first 
moment, (E^ = | dE i E i P q {E), defined in terms of the escort distribution [9] associated with 

m 



P(E) = (41) 

P ' {E ' ] - r^mr 



is finite if q is chosen to be 



4 = 1 + -. (42) 
s 



This is the minimal requirement for the finiteness of (E ( ) . The denominator in eq. (41) and 
(is,-) are calculated to be 

c,=]dE, [/(£,)]' J^-^{fi), (43) 



».-<*>. = 7*?- (44) 



respectively. Thus, the constant a in eq. (36) [and eq. (38)] is given in terms of the physical 
quantity u q , i.e., the generalized internal energy. From eqs. (43) and (44), we have 



c q =—u\-\ (45) 



To connect these quantities with those of thermodynamics, (5 and u q are regarded as a conjugate 

pair in the Legendre transform structure [10]. For this purpose, an entropy functional, S q [f], is 
defined such that 
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(47) 



where 



r< G Wf -£ V (48) 



The superscripts "(G)" are used to indicate the quantities derived within the scope of the 
generalized central limit theorem. The relation between £ and (5 appearing in eq. (34) can be 
established, once the form of the entropy S^ G) [/] is specified. At this stage, we examine the 

following form for the entropy: 

Ct/H^K^- 1 )' (49) 

where a scale o is introduced to make S^ G) [/] dimensionless. Using eq. (46), we find 
v u~ q 

P = — ?-. (50) 

From eqs. (48) and (50), we ascertain the validity of eq. (47). 

It may be mentioned that in recent papers [3,4] similar studies on convergence properties in 
both the full and half spaces have been studied. 



Principle of Maximum Tsallis Entropy 
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Actually the entropies appearing in eqs. (25) and (49) and the distributions in eqs. (16), (30), 
and (34) are the Tsallis nonextensive entropy [13] and its optimal distribution under the 
constraint on the generalized first moment with respect to the escort distribution in eqs. (18) (41), 
respectively. 

To complete our program, here we make contact of the above three derivations with 
nonextensive statistical mechanics. This is done by maximizing the Tsallis entropy [13] 

^W= T ^jEW«)]*- 1 j < 51 ) 

1 I m J 

subject to the constraints on the normalization of the probability p(m) and on the normalized q- 
expectation value of a physical quantity, a(m) [14]: 

£ a(m)[p(m)Y 

(A) = " ^— = a. (52) 



Introducing the Lagrange multiplier (3 to take into account the constraint in eq. (52), we obtain 
the following results: 



ZTW ^ ^ (53) 



Z\ M \P) = Y j e q {-P[a{m)-a\ 



l-q 



(54) 



(55) 



(56) 
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Here, the superscripts "(M)" is used to indicate the quantities obtained from the principle of 
maximum Tsallis entropy. Note that 



q yH '=0, (57) 



in view of eq. (52). We now define T^ M \f5) by means of the Legendre transform 

r^(/3) = \n q (zM(p))-Pa, (58) 



so that 



dp 



= —a. 



(59) 



It is straightforward to show that 

W[p,} 



da 



P- (60) 



Note that eqs. (59) and (60) define a canonical pair, which is of central importance in 
thermodynamics. It is useful to note that the above results go over to the standard ones in 
Gibbsin canonical ensemble theory in the limit q — > 1. 

In conclusion, we have shown that the power-law canonical distribution can be derived within 
the fundamental principles of statistical mechanics based on four different methods: steepest 
descent method, counting method, the generalized central limit theorem, and the principle of 
maximum Tsallis entropy. 
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